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Preliminary results from a Fermilab study of the sensitivity for higgs production at the Tevatron 
in run II and beyond are presented. The study extends existing results by systematically combining 
results for all decay channels, considering the production of higher mass higgs bosons and interpreting 
the results in the context of supersymmetric higgs production as well as standard model production. 
In addition new analysis methods which significantly improve sensitivity are used. 

ON '. I. INTRODUCTION 

ON ■ 
ON 1 

— , ' The standard model of particle physics has been studied with very high precision over the course of the past ten years, 
£>^' and no significant deviations have been found. Despite this, our understanding of the origin of electroweak symmetry 
breaking is still incomplete. This arises in large part because the only remaining undetected standard model particle, 
| the higgs boson, mediates electroweak symmetry breaking in the standard model. The highest available center-of- 
qq ' mass energy for the years 2000 to 2004 will be at the Fermilab Tevatron pp collider with y/s = 2.0 TeV. It is natural to 
explore the sensitivity to higgs production at the Tevatron. This paper contains preliminary results from a year-long 
, study conducted jointly by the Fermilab theory group and the CDF and D0 experiments.^] The goal is to quantify 
^ ' the higgs discovery potential at the Tevatron in the coming run II and possible extensions. Results are presented as 
the luminosities required to exclude higgs at the 95% confidence level, or to establish either 3er or 5er excesses over 

■ predicted backgrounds. 

The starting points for this study are the higgs mass constraints expected from LEP2 H and previous Fermilab 
studies H ||. This study extends the previous Fermilab results by (1) including additional standard model decays in 
ON ■ the mass regions previously explored, (2) testing the sensitivity for higgs masses M H > 135 GeV, (3) systematically 
combining results from all channels, (4) interpreting the results as supersymmetric (SUSY) higgs production and (5) 
Qh! considering additional decay modes arising from SUSY models. In addition, a detector simulation was developed 
O ,' which gives significantly more realistic event reconstruction than some of the previous studies used. 

This paper has six sections. The first describes the production and decay of standard model higgs bosons and the 
• • , simulations used in this study. The second and third sections contain results for standard model higgs production 
.£h , m the mass ranges 90 < Mh < 135 GeV and 135 < M# < 200 GeV respectively. The fourth section presents the 
combination of the results in sections two and three. The fifth section describes the extension of the results to SUSY 

■ production. The last section describes studies of additional SUSY-specific decays, particularly final states having four 
&-quarks. 

II. PRODUCTION, DECAY, EVENT GENERATION AND DETECTOR SIMULATION 

The production cross sections and decay branching ratios for a higgs bosons have been calculated by a number of 
groups. H Those for a standard model higgs boson are shown in Fig. [j]. These plots indicate that the highest cross 
section production modes are pp — > H , pp — > WH and pp — > ZH . The higgs decays dominantly to the most massive 
kinematically allowed final state. For Mjj < 135 GeV, the dominant decay mode is H — > bb with a branching ratio 
of roughly 80%. For Mh > 135 GeV, the dominant mode is H — » WW. Thus, searches for lower-mass higgs will 
be looking for final states with at least two ^-flavored jets, and the higher mass searches will have multiple (virtual) 



x For more information, see http://fnth37.fnal.gov/susy.html 
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W bosons. For most of the mass range in question, the pp — > H mode has very poor signal-to-noise, and the most 
useful modes are the pp — > WH and pp — > modes with the W ox Z decaying to leptons. 
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FIG. 1. Production cross sections and decay branching ratios for the standard model higgs boson as a function of mass. 

Unless explicitly noted, events used in these analyses were generated using the Pythia ||, Isajet || or CompHep 
programs. The generated four-vectors were then input to a detector simulation program, SHW, developed for the 
run II workshop § . This program uses parameterized resolutions for tracking and calorimeter systems and particle 
identification to perform simple reconstruction of tracks, calorimeter-based jets, vertices and trigger objects. The 
resolutions used represent a typical run II detector and are drawn from CDF and D0 internal studies. Particle 
identification efficiencies are included by parameterizing results from other CDF and D0 studies. 

The SHW program was verified by comparing selection efficiencies between SHW and data or between SHW and 
well-established run I simulations used by CDF or D0. The most stringent test was a comparison of nearly identical 
analyses of the pp — > WH — > (lv)(bb) channel. Two analyses of this channel have been performed, one based on 
a run I CDF simulation with the geometrical acceptance extended to correspond to the run II CDF detector and 
the second based purely on SHW. The first analysis predicts 5.0 signal events and 62.8 background events/fb^ 1 for 
Mh = 110 GeV. The second predicts 4.5 signal and 62.5 background events for the same conditions. 

III. LOW MASS HIGGS SEARCHES, M H < 135 GEV 

When Mh < 135 GeV, the dominant decay mode is H — ► 66. Analyses have been performed for all pp — ■> WH 
and pp — > ZH final states.^ The possible final states are: (1) pp — > WH — * ivbb, (2) pp — * ZH — * vvbb, (3) 
pp — > ZH — > £ + £~bb and (4) pp — > WH — > qqbb or pp — * ZH — * qqbb. The primary backgrounds to these channels are 
W + bb and Z + 66 with the 66 pair from gluon radiation, single top-quark production and top-quark pair production. 

All analyses for these channels begin with a preliminary selection based on the number and type of final state objects. 
For example, the pp — > WH — > tvbb analysis requires a charged lepton with Et > 20 GeV, missing transverse energy 
fr > 20 GeV and two 6-tagged jets having Et > 15 GeV. Similar selections are applied for the other channels. After 
the basic selection, a requirement is made that the mass of the reconstructed 66 system be within (typically) 2a of 



2 The mode pp — » H — » bb was considered, but the signal to noise was too poor for it to have any sensitivity when compared 
to the pp —> WH and pp — > ZH modes. 
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the generated higgs mass. Additional clean up requirements are also made. As an example, in the pp — ■> ZH — > vvbb 
channel, there can be no isolated tracks with px > 15 GeV. This rejects events with high-pr leptons which failed the 
lepton identification. The resulting number of signal and background events corresponding to 1 fb _1 of data are given 
in table [| The pp — > WH — ► ivbb and pp — ► Zi7 — > modes offer the best sensitivity with the pp — > — > £ + £~bb 
mode not far behind. The all-hadronic final state looks quite difficult. 

In addition to these analyses, a multivariate analysis using neural networks has been peformed for the pp — ► WH — ► 
Ivbb channel. This style of analysis has been used with considerable success by D0 in the top mass Q and all-hadronic 
top decay analyses. The basic principle is to exploit correlations within an event in an automatic manner. The 
left panel of Fig. |^ shows the number of predicted signal and background events for pp — > WH — ► ivbb analyses. 
Each point in the figure represents one possible analysis. The band labelled "rgsearch" corresponds to hypothetical 
analyses performed using selections using the standard technique of sequential requirements applied to event variables, 
with each requirement a single-valued comparison such as fr > 20 GeV. The point labelled "TeV 2000" is the result 
from a previous Fermilab study 0. The point labelled "neural net" is the result from the multivariate analysis. One 
sees that for a fixed background, the signal is increased by roughly 50% using the neural network. Similar gains are 
expected in all other channels in this mass range. 
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FIG. 2. The left panel shows predicted signal and background selection efficiencies from the neural network analysis. Each 
point corresponds to a possible event selection. The point labelled "neural net" is the result from the multivariate analysis 
described in the text. The right panel shows neutral network output equal probability contours in the Ht vs. Mjj plane. Ht 
is the scalar sum of all jet energies, and Mjj is the invarient mass of the tagged dijet system. The open boxes are background 
events, and the closed boxes are signal. 



A. Other Improvements 

Results have also been obtained for hypothetical improvements in mass reconstruction and 6— jet tagging efficiency. 
The analyses were repeated after artifically improving the reconstructed dijet mass resolution in steps up to a 50% 
better resolution. The results in Table | include an improvement in mass resolution of 30%. This level of improvement 
is possible when information such as charged track energy is used in the mass reconstruction in assocation with the 
calorimeter-based jet energies currently used. Such an improvement has already been realized in a preliminary run I 
CDF analysis of the pp — > Z — > bb channel. |ll[] Improved mass resolution offers considerable benefits because for a 
selection with a fixed signal expectation, the background will decrease as the resolution improves. 
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The effect of improved 6— jet tagging has also been explored by artificially improving the second jet tagging efficiency 
by up to a factor of two. The gains from this improvement are not as important as those from mass resolution 
improvements because both signal and background increase with improved tagging efficiency. 

IV. HIGH MASS HIGGS SEARCHES, M H > 135 GEV 

Previous Fermilab studies have concentrated on the lower mass higgs states which decay dominantly to bb. This 
study includes analyses designed for final states in which the higgs decays to WW or ZZ instead of bb. This corresponds 
approximately to Mh > 135 GeV. Three final states are considered: (1) Three leptons, l^l'^P, arising primarily 
from pp — > WH — > WWW, (2) Dilcptons and neutrinos, l + l~vv, from pp — > H — > WW and (3) Like-sign dilcptons 
plus jets, l^l^jj, from pp — > WH — ► WWW and pp — ► Zff — * ZWW. |E| The dominant backgrounds are standard 
model production of WW, WZ, ZZ, and WiZ) + jets and it and multijet events with misidentification arising from 
detector effects. The standard model sources dominate the detector effects.^] 

As for the low mass analyses, the initial selections are based on simple variables related to the boson decay-product 
kinematics. However, to reach usable sensitivity, the analyses then use either (1) requirements typically relating to 
angular correlations arising from spin differences between signal and background or (b) likelihood methods. In both 
cases new variables have been designed. Figure |^ shows one such variable used in the l + l~vi> analysis, the cluster 
mass Mq = yP^Jt) + m 2 (££) + A result of the tuning is that the signal and background have similar mass 

distributions, so these analyses must be treated as straight counting experiments. The numbers of expected signal 
and background events for the high-mass channels are given in table ||. 
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FIG. 3. The cluster mass variable for background (shaded region) and for signal and background together (open region) for 
the l + l~uv analysis. 

V. COMBINATION OF STANDARD MODEL SEARCH CHANNELS 

The results in the preceeding two sections have also been combined to form a single unified result. Figure |i| shows 
the luminosities required for 95% CL exclusion, 3<r evidence and 5a discovery as a function of standard model higgs 



3 In general, the backgrounds arising from detector effects use conservative misidentification probabilities based on run I 
analyses by both experiments. 
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mass. These contours include statistical and systematic errors^] and the channels are combined using the prescription 
of reference Ol . 
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FIG. 4. Luminosity required to achieve 95% confidence level exclusion, 3cr evidence and 5cr discovery as a function of higgs 
mass. The results use all channels and assume results from both CDF and D0 having equal sensitivity. The experimental 
uncertainties used are described in the text. 



VI. SEARCH FOR NEUTRAL HIGGS BOSONS IN SUSY MODELS 

Most supersymmetric extensions of the standard model result in five physical higgs states, two neutral scalars 
denoted h and i?(with Mh < Mh), a single pseudoscalar A and the charged doublet H ± . Unlike most other 
supersymmetric particles, the masses of the higgs particles depend only on two parameters, typically chosen to be 
tan/3 and Ma (or M/,). Furthermore, the h boson must satisfy Mh < 130 GeV, and it may easily be less massive 
that this. For much of SUSY phase space, the h decay is identical to that of the standard model higgs, and only the 
production cross section differs. Given this, the searches for low mass standard model higgs can be easily converted 
to searches for SUSY higgs. The left-hand panel of Fig. || shows the 5a exclusion contour in the tan/3 vs. Ma plane 
for the case in which all non-higgs SUSY masses are around 1 TeV, the systematic error is 10%, and CDF and D0 
data are combined. One sees that for this example, SUSY models having 1 < tan (3 < 50 and 80 < Ma < 400 GeV 
are excluded. 



4 The systematic errors are assumed to scale with luminosity. The scaling is expected to hold at least until 2% relative 
systematic errors are reached. Systematic uncertainties at this level do not limit the analyses. 
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FIG. 5. Exclusion contours for SUSY models. For the left panel the shaded region is excluded based on a reinterpretation 
of the standard model searches in section two as SUSY higgs production. This example shows the 5a exclusion region with an 
integrated luminosity of 20 fb -1 . The right panel shows the result from the pp — > <f>bb — > bbbb analysis described in section VII. 

VII. SEARCHES FOR THE SUSY MODE pp -» </>&& -> bbbb, <j> = H, A, H 

Within SUSY models, new higgs production modes exist, some with couplings proportional to tan 2 (3. For sufficiently 
large values of tan/3 these become the dominant production modes. One such mode is pp — > hbb^. This results in 
final states containing four 6-flavored jets. Analyses of this channel have been carried out. These typically require 
four jets, three of which satisfy &-tag requirements. All possible mass combinations of 6-jets are computed, and the 
resulting distribution is examined for a peak near the generated higgs mass. The 95% C.L. exclusion contours in the 
tan/3 vs. Ma plane are shown in the right panel of Fig. ^ for a variety of integrated luminosities. 

VIII. CONCLUSIONS 

Studies of the experimental sensitivity to higgs production for Tevatron Run II and beyond have been carried out. 
Both standard model and supersymmetric higgs production have been considered. It is found that with 4 fb -1 of data, 
standard model higgs can be excluded at 95% confidence over the interval Mjj < 125 GeV and 155 < Mb < 175 GeV. 
With 10 fb -1 , a standard model higgs boson will be seen as at least a 3<r excess over the mass ranges Mh < 125 GeV 
and 145 < Mh < 175 GeV. These results have been converted to limits on the SUSY parameter space, with one 
example shown. 
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TABLE I. Numbers of expected signal and background events for each low-mass channel in 1 fb .A 30% improvement in 
mass resolution over that from SHW has been assumed. See the text for details. 
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Channel 


Rate 


90 


100 


110 


120 


130 




S 


2.5 


2.2 


1.9 


1.2 


0.6 


vvbb 


B 


10 


9.3 


8.0 


6.5 


4.8 




s/Vb 


0.8 


0.7 


0.7 


0.5 


0.3 




S 


8.4 


6.6 


5.0 


3.7 
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Ivbb 


B 


48 


52 


48 


49 


42 




s/Vb 


1.2 
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0.5 


0.3 




S 
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0.5 


0.3 


l+rbb 
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3.6 


3.1 
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1.1 
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0.5 


0.5 


0.4 


0.3 
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8.1 


5.6 


3.5 
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6800 


3600 


2800 


2300 


2000 




s/Vb 


0.10 


0.09 


0.07 


0.05 


0.03 



TABLE II. Numbers of expected signal and background events in 1 fb 1 for the high-mass channels. 
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l^l'+P 
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0.73 


0.73 


0.73 


0.73 


0.73 




s/Vb 


0.05 


0.09 


0.13 


0.14 


0.18 


0.12 


0.11 
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2.6 


2.8 
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1.1 


1.0 


l + Tvv 


B 






44 


30 


4.4 


2.4 


3.8 




s/Vb 






0.39 


0.51 


0.71 


0.71 


1.9 




S 


0.1 


0.20 


0.34 


0.53 


0.45 


0.38 


0.29 




B 


0.85 


0.85 


0.85 


0.85 


0.85 


0.85 


0.85 




s/Vb 


0.11 


0.22 


0.37 


0.57 


0.49 


0.41 


0.31 
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